Motion of Colored Particle in a 
Chromomagnetic Field 

Sh.Mamedov* 

High Energy Physics Lab., Baku State University, 
Z.Khalilov str.23, Baku 370148, Azerbaijan 
and 

Institute for Studies in Theoretical Physics and Mathematics (IPM), 
P.O.Box 19395-5531, Tehran, Iran 



Abstract 

The Dirac equation in a chromomagnetic field is solved for colored particle moving 
in a limited space volume. Quantized energy levels and the corresponding wave func- 
tions are found for backgrounds both directed along third axes and having spherical 
symmetry. It was shown interrelation with the case of motion in an infinite space 
volume. 
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Introduction 



Quantum and classical mechanics of non-abelian charged particles has been studied by var- 
ious authors. Problems in this theory serve to find the solution of QCD problems and are 
connected with description of quark's motion inside hadrons or any other limited space vol- 
ume. Since QCD has confinement property, presumably it will be better to consider QCD 
problems in a finite volume as well [1]. As is known, there are color background fields in 
the QCD vacuum [2] and so, it makes sense to study the problem of motion of a colored 
particle in a color background in limited space volume. The simplest case, a constant color 
background could be given by two different types of vector potentials [3,4]. Problems in 
an external color field given by abelian vector potentials are solved analoguosly to the ones 
in abelian theory. The second type of vector potentials are non-commuting constant vector 
potentials, which are not gauge equivalent to first type [4] . These potentials are used to solve 
the ground state problem of QCD [4,5] and other different problems connected with the QCD 
vacuum [5-9] as well and sometimes give physically different results from that ones in which 
first type potensials was used. Here we shall solve Dirac equation in a finite space volume 
in a constant background chromomagnetic field given by the second type vector potentials. 
The solutions and spectra found are more suited to the physical situation and by use of 
them could be constructed quark's Green's function by the help of exact solution method 
[10], which could be used in solving QCD problems. 



1 Constant homogenous background 

Let us define an external chromomagnetic field by constant vector potentials. Within the 
SU(3) color symmetry group they look like 

A\ = yfr8 la , A\ = ^2a, A a 3 = 0, Al = 0, (1.1) 

where a = 1,8 is a color index, r is a constant and 5^ a is the Kroneker symbol. The field 
(1.1) is directed along the third axes of ordinary and color spaces : 

F* 2 = Hl = gr other F« u = 0. (1.2) 

Here g is the color interaction constant. 

The Dirac equation for a colored particle in the external color field has the form 

( 7 ^P M -M)^ = 0, (1.3) 

where P M = p^ + gA^ = p^ + gA®^- , A a are Gell-Mann matrices describing particle's color 
spin. Equation (1.3) written for Maiorana spinors <p and x nas t ne well-known form 

(^)V = "(^ + M 2 ^, (1.4) 

where the Pauli matrices a 1 describe a particle's spin. Here and afterwards -0 means or x- 
The spinors <p and x have two components, corresponding to the two spin states of a par- 
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tide ip = 



. Each component of ip transforms under the fundamental representation 



of color group SU(3) and has three color components describing color states of a particle 
corresponding to three eigenvalues of color spin A 3 



± 



<MA 3 = 




( v4 1} \ 








<MA 3 = 


o) J 


{ 4 3) J 



1.5) 



Since the background field (1.1) is time independent, the equation (1.4) has got the form 
T + 9 ^h + gr* (pxA 1 +p 2 A 2 =F ^A 3 ) V± = {E 2 - M 2 ) (1.6) 



Here I 2 = 



1 

1 ] is the color matrix. Writing this equation for the color components 


'1.5) we get two independent systems of differential equations 

(pi - ip 2 ) i> {2) = {E 2 - M 2 ) ^ 
J? 2 + g 2 r) Vi 2) + gr^ (pi + ^p 2 ) = {E 2 - M 2 ) ^ 

[E 2 - M 2 ) ^ 3) 



-^1 /(3) 



(p> 2 + g 2 r) 4> { y + gr\ ( Pl - ip 2 ) ^ ] = (E 2 - M 2 ) 



(i) 



:i-7) 



-ft 2 iP {2) + ^ (pi + ip 2 ) 4> { - ] = (E 2 - M 2 ) ^ 
-f 2 ^ = (E 2 - M 2 ) ^ 

From (1.7) we get the same equation for all states ip^ (i = 1,2) 



J? 2 + 9 -^ + M 2 -E 2 ) - 9 2 r(p 2 ± + 9 -^ 



0, 



(1.8) 



where p\ = p\+p 2 - This equation possesses rotational invariance around the z axis. If we 
consider free motion of the particle in all a? space, then we should not impose any boundary 
condition on solution of (1.8) and shall get continous energy spectra found in [11,12] for 
plane wave solutions ip^ ~ e l P r 

2 

(1.9) 




g 2 r gr-< 



+ p 2 + M 2 . 



This spectrum is used in solving of various problems [5-9]. Since free motion of the colored 
particle in nature could be take place in limited space, for example inside a hadron, let us 
consider motion of particle in limited space bounded by a cylinder with radius r and hight 
z and solve (1.8) on boundary conditions tp± (r = r , z) = 0, ip± (r, z ) = (r 2 = x 2 + y 2 ) . 
Let us rewrite (1.8) in the form 

, 2 



(-V 2 + 4 I ) -2(E 2 -M 2 + 
(E 2 -M 2 -^f -g 2 r 



9L 

dz 2 



/(») 

V4 



-V 2 + ^ 
0. 



)+ 



(1.10) 
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With (— V 2 + ^j-^j tp± = r(± it has got the following form 

(-v 2 + £f) r,f -2(e 2 - m 2 + £i) n f+ 



Acting on equation (1.11) by operator (— V 2 + we get an equation for r)± that has the 

same form as for ip± , i.e., (1.10) with the replacement ip^ on r)±. This means that 77^ and 
ip^ differs only by a general constant multiplier k 2 : 77^ = k 2 tp± or 

(-v 2 + ^)^? = fcV£°- (L12) 

Equation (1.12) is equivalent to (1.10). In the appendix is shown another way to reduce 
(1.10) to the form (1.12). In a cylindrical coordinate system the solution of (1.12) can be 
looked for using the separation ansatz ip± ("f > ) = ip (r) • u (<p) ■ v (z) and (1.12) is divided 
into three independent equations 

-*0± = m 2 u^) (1.13) 
^ + AM,)=0 

The solutions of last two equations with the boundary condition v (z ) = are u ((f) = 
-j=e im v (m = 0,±1,±2, ... and v (z) = sinp 3 z (A 2 = p 2 , p 3 = ^) . For the p 3 = case 
we can choose v (z) = 1. It is worth to remark that (1.12) and (1.10) become equivalent 
after taking the last equation in them into account and there is no such problem for the 
case p 3 = 0. Taking into account explicit expression of v (z) and (1.12) in (1.10) we get the 
following equation for the constant k 2 



,2^\ / Jl 



2 



(k 2 ) - 2k 2 [E 2 - M 2 + + [e 2 - M 2 - + g 2 rp 2 3 = 0, 

from which one finds the following expression for the constant k 2 
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*? >2 = \^E 2 -M*-p 2 ±^j +pl (1.14) 

The first equation in (1.13) is the Bessel equation. Under the finitness condition on r — > 
it has a solution determined by the Bessel function 



Mr) = J m ^kl 2 -p 2 - 9 ^ . (1.15) 

If we do not impose any boundary condition on the states (1.15) with given m, then they 
have an energy from continous spectrum (1.9), which does not depend on m. Imposing the 
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boundary condition ^i j2 (r ) = we find two branchs of quantized energy levels of the colored 
particle in a chromomagnetic field (1.1) 



\ 



v 2 

(N)\ 2 _ i\ 



o 2 r ot2 



+ Ps + M J . (1.16) 



/ 



Here affl are the zeros of Bessel function and the quantity N labels the sequence of zeros 
N = 1,2,3.... From (1.16) one sees that these energy branchs, as branches of the continous 
spectrum (1.9), are defined by ± sign of the field strength, which is not connected with 
direction of field. The energy levels E^p are determined by the m -chromomagnetic quantum 
number, i.e. by the projection of the chromomagmetic moment of the particle onto the 
chromomagnetic field. Another quantum number is N, determining the energy level 
a result of quantization due to the finitness of the motion. In contrast to the discrete 
spectum (1.16) the continous spectum (1.9) is not determined by the quantum number m. 
Thus, finitness of space volume turns continous spectra (1.9) for a given m into a discrete 
series determined by the quantum number N. Both spectra (1.9) and (1.16) are infinitely 
degenerate. Since (1.12) is universal for the spin ± and color (i) indices, any of these states 
may have any energy from spectra (1.16) or (1.9). If in the energy branchs (1.16) one makes 

(JV) 

the replacement = p± we get continous spectra (1.9). This is suitable to quantization 
of momentum in standing waves. Taking into account spectrum (1.16) in (1.14) we find the 
value of constant k 2 , 



which is the same for both energy branchs and so, the solution (1.15) is the same for these 
branchs. Thus, the wave functions ip^ ("r*') from chosen energy branch ^Effl^j or (jZffl^j 
are equal to 

^ +°° 1 / a (iV) \ 

1>£(?) = ]T e ^sm(p 3 z)J m r^- . (1.17) 

m=-oo V 27T V r / 



Thus, the motion of colored particle in the field (1.1) along z axis is a standing wave and 
in the (x, y) plane is a rotation. Writing (1.13) for the radius a of the cirle in which particle 

rotates 

1 d ( dip{a)\ ( 2 g 2 T 2 m 2 \ 
ad-r{ a —) + [ k - — -P3-^)^^ = ° 

we can find discrete values for these orbits 

(N) 171 

a( m } = —ro- 
am 

In the continous spectrum case the constant k 2 and radius a are 

k = p + — , a = — 
4 p ± 
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and the solution (1.15) becomes 



^1,2 (r) =ip(r) = J m (rp±) ■ 



The states ip^ (T*) correspond to states of freely moving colorless scalar particle. 



2 Spherically symmetric configuration of background 

Let us consider the case of spherical symmetry in the ordinary space with constant chromo- 
magmetic field defined by non-commuting potentials 

A\ = V^Sia, A a 2 = ^8 2a , Al = v^ 3 a, A a = 0. 
The non-zero components of the field strength tensor are 

rpl p2 1713 

^23 ~~ r 31 ~~ r 12 — y' ■ 

In the external field (2.1) the explicit form of (1.4) is 



~]? 2 + M 2 + + grh a p a - ^-a a \ a \ # = E 2 ^ . 



(2.1) 



(2.2) 



Let us denote components of the spinor \l> by ipi and ip 2 '■ * 



4>2 



. Since the external 



field has three components in both spaces, in this case, different from the previous one, 
the components ijjfl do not describe states with definite value of the spin and color spin 
projection onto the field. Equation (2.2) can be written in terms of a system of differential 
equations for the components ip^l 

(A + grhs ~ 4 1 ) *Pi ] + 9^ ( P i - ip 2 ) 4? = 
[A - gAp 3 + ^) 4? + gA (pi + ip 2 ) 4? = 5 2 T'02 1) 

(t 2 + M 2 ) 4? = E 2 4 3) 
[A + gAp 3 + ^f) 4^ + gA (pi - ip 2 ) ip { 2 2) = g 2 T%l) {2) 
(A - gAp 3 ~ 4 1 ) 4 2) + 9t\ ( Pl + ip 2 ) 4 2 1] = 
(J 2 + M 2 ) 4 3) = E 2 4 3) 

where operator A denotes A = ~p" 2 + M 2 + — E 2 . From the system (2.3) we get the same 



(2.3) 



differential equation for the all states 4 2 (i = 1,2) 



4% 



0, 



(2.4) 



which possesses rotational invariance. Since the operator in the first square brackets com- 
mutes with the second one, (2.4) could be divided into two equations 

,2 

2 ' 

2 



2 — >2 

g tp 



*8 



0, 



(A + &) -g 2 r(t 2 + 9 2 r) 



0. 



(2.5) 
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From (2.5) we get four branchs of the continous energy spectrum found in [11] for plane 
wave solutions ip±\ ~ e* ^ r 



E 2 

^1,2 



1 \ 2 
r3 



(2.6) 



In the same manner as used in the previous section from (2.5) we get the equivalent equations 

(-V 2 ) # = tf^f, (" V 2 +A) V>f = ftf^j = 1,2). (2.7) 

If we take (2.7) into account corresponding equations (2.5) we found values of k 2 d 



k\ 2 = k\ A = k 2 = VE 2 - M 2 ± 



gr 7 < 



(2.8) 



The solution of (2.7) can be found in textbooks of quantum mechanics ( [13] ). In a spherical 
coordinate system for the wave functions tp^ one used separation ansatz, i.e. it is chosen as 
a multiplication of the radial R(r) and angular Y™ (9, (p) parts 

# m = R(r) ■ YT (0, <p) ■ 



Here r = y/x 2 + y 2 + z 2 , / and m are the orbital angular momentum and chromomagnetic 
quantum numbers, 9, ip are the polar and azimuthal angles, respectively. The spherical 
functions Y™ (9, ip) are expressed by means of Legendre polynomials Pj"^ (cos 9) 



Yl n (9,p) = a r , 



(21 + 1) 


(l- 


m \)\ 


47T 


(1+ 


m \)\ 



PP (cos9)e 



%rmp 



(2.9) 



with 

Qj Tl 



1 on m > 



m < 



The equations for the radial parts of wave functions are the same as for many quantum 
mechanical problems possessing rotational invariance 



V 2 r R(r) + [kl- 



1(1 + 1) 



R(r) = 0, 



(2.10) 



where v = 1,2,3,4; kf 2 = k 2 , k\ A = k 2 - g 2 r. With the notations Q (r) = ^frR(r) (2.10) 
turns into Bessel's equations for Q (r) 



Q" (r) + -Qf (r) + 

r 



1C - ^ + }' \ O(r) = 0. 



(2.11) 
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The function R(r) must be finite on r — > 0. This means that solution of (2.11) are first type 
Bessel functions 

RZ = -£^J l+1/2 (k v r). (2.12) 

Imposing boundary conditions Ri{vq) = we enclose motion of the particle by a sphere 
with a radius r , where r agrees with the hadron's size. Since, the external field (2.1) 
is not depend on r, we have got the same equation for R(r) as for a freely moving particle 
enclosed in a sphere, differing only from the expression by k 2 constant. From these boundary 
conditions one finds quantized energy levels of spectrum branchs 



(2.13) 



In this case the boundary condition on R\(r) with fixed value of the angular momentum 
quantum number I chooses from the continous energy spectrum (2.6) the discrete series 
(2.13) of the radial quantum number N for this value of /, and the energy degeneracy for 
this state remains infinite. The replacement connecting discrete and continous spectra is 

(N) 

=| J? |. By means of spectra (2.13) one finds the same values for k u 

(AT) 

k ~ k W ~ ^ 

The radius a of the particle could be found from (2.10) using the maximum condition R'(a) = 
and is equal to 

AN) _ r Vj (j + 1) 

As the energy spectrum is described by (2.13), the motion orbits of the particle are quantized 
and defined by the orbital and radial quantum numbers I and iV as well. Thus, general 
solution of (2.7) is 

oo I 

^ = E E CTYT (0, <p) -jaf-J, + i /2 (kr) ■ (2.14) 
i= 0m =-i a\ V r 

In the case of this field configuration the particle moves on "s", "p","cT, "/".... orbitals 
known from orbital angular momentum eigenfuctions. For the continous spectrum the con- 
stant k and radius a are 

I P I 

The plane wave could be decomposed by means of spherical waves and therefore, plane wave 
solution obey initial equation as well [13]. 

Using the energy spectra (1.16) and (2.13), estimates for field strength of gluon condensate 
[2] and hadron's size the energy of a photon emitted by excited hadron may be calculated. 
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Using the wave functions (1.17) and (2.14) the explicit expressions of Green's functions for 
these particles in finite temperature can be found [10], which gives a possibility to investigate 
different radiation corrections in the chromomagnetic field and the quark condensate 



< ^ >= -i lim Tr [G (z) - G (z) \ A =o\ ■ 



z->0 



(2.15) 



The spectra found, (1.16) and (2.13) and the quark condensate (2.15) open new possibilties 
to investigate chromomagnetic catalysis of color superconductivity as well [6] . 

Appendix 



It is easy to reduce (1.10) to the form (1.12) in momentum representation 



rl>g (?)=[' —< 
J-oo (2tt) 



^ (?) . 



Then (1.10) will have the form 



2 / 2 

2 t?2 \ ,2 1,2 , 9 T 



1*1' (t) = 0- 



which is divided into two equations 



(?) = o 



and can be written in the following form 



(A.l) 



Ipi + 9 -{ ± ^ ] m = - m 2 -pi) ® (?) . 



(A.2) 



From (A.2) we get 



1/2 



Squaring the last equality and returning to the coordinate representation by an inverse 
Fourier transformation we get equation (1.12) 



e 2 -m 2 - p It — } 4>±- 
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